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We present the analysis of the distribution of the distance of muon pairs detected deep underground by the 
MACRO experiment. A detailed analysis of the low separation region pointed out an excess of experimental data 
with respect to the Monte Carlo expectation. The effect is discussed taking into account the /^+N— >/u+N+/i + /i~ 
process. 



1. Introduction 

The distribution of the muon pair separation 
deep underground (decoherence function) plays 
an important role in the study of Cosmic Ray 
physics in underground experiments, being sen- 
sitive to some aspect of the hadronic interac- 
tion model. Taking advantage of the large area, 
MACRO 0] can access the largest separation, ex- 
ploring thus the highest P t regions. Moreover, the 
capability of the MACRO detector to resolve very 
close tracks permits the extension of the decoher- 
ence analysis to a distance region hardly studied 
in the past. 

2. Analysis of the behaviour at small dis- 
tances 

We have analyzed about 3.4 ■ 10 5 events, corre- 
sponding to 7732 hr live time using the lower part 
of the apparatus. These events were submitted to 
the following selection criteria: 
- Zenith angle smaller than 60°: this choice is 
dictated by our limited knowledge of the Gran 
Sasso topographical map for high zenith an- 
gles. Furthermore, we cannot disregard the at- 
mosphere curvature for larger zenith angles, while 
at present our simulation model does not include 



*For the complete list of the Collaboration see the paper 
" Relevance of the hadronic interaction model in the inter- 
pretation of multiple muon data as detected with MACRO 
experiment" by O. Palamara at these proceedings 



that. 

- The number of streamer tube hits out of track 
smaller than 45. This selection aims to eliminate 
possible misleading track reconstruction in events 
produced by noise in the streamer tube system 
and/or e.m. interactions inside or in the neigh- 
bourhood of the apparatus. 

- Track pairs which survive a loose parallelism 
cut: this rejects hadrons from photonuclear inter- 
actions close to the detector. This selection also 
eliminates tracks reconstructed using hits due to 
e.m. interactions in events survived to the previ- 
ous cut. 

The last cut is not completely efficient in re- 
jecting muon tracks originating from local parti- 
cle production because the angle between these 
tracks may fall within the limits imposed by the 
parallelism cut. Moreover, these limits cannot be 
further reduced since the average angular diver- 
gence due to multiple muon scattering in the rock 
overburden, is about 1° at the MACRO depth. 

A further selection was applied in order to re- 
duce these effects. We computed, for each muon 
track in the wire view, the ratio R between the 
number of streamer tube planes hred by the muon 
to the number of expected fired planes consider- 
ing the track direction. Only tracks with R > 
0.75 were accepted. The application of this cut 
(hereafter the C4 cut) in the wire view alone is a 
good compromise between the rejection capabil- 
ity of the algorithm and the loss of events due to 



3 

< 1 

S -1 

z 10 

■D 

9 -2 

- 10 



10 



10 



10 



10 



10 



• DATA 

O HEMAS (MACRO fit) 



\ 



V 



% 



h 



'ff + - 



\\ 



20 



40 60 

Separation (m) 



Figure 1. Comparison between the experimental 
decoherence function with the HEMAS expecta- 
tion. 

the unavoidable inefficiency of the streamer tube 
system. We found that in the wire projective view 
the probability to reject a muon track for a casual 
alignment contiguous inefficient planes is 2.0%. 

To test the ability of this cut to reject hadronic 
tracks, we used FLUKA || to simulate 3028 hr 
of live time in which muons were accompanied by 
hadronic products of photonuclcar interactions in 
the 10 m of rock surrounding the detector. We 
found that the parallelism cut alone provides a 
rejection efficiency of about 54.6% of the pair 
sample, while the addition of cut C4 enhances 
the rejection to 95.9%. Furthermore the effect of 
hadron contamination, is very small, contributing 
less than 1% in the overall muon pair sample. Af- 
ter the overall application of these cuts, the num- 
ber of unambiguously associated muon pair tracks 
survived is 355,795. Fig. 1 shows the compari- 
son between real data and HEMAS expectation: 
a general good agreement is found, showing the 
HEMAS capability in reproducing the muon bun- 
dle features deep underground. Nevertheless a 
relevant discrepancy at low distances is observed, 
indicating a possible physical source for this ex- 
cess in the real data. For istance, considering 
separation smaller than 80 cm, an excess of 55% 
is observed in the real data(Table 1). The pro- 



cess of muon pair production by muons in the 
rock is a natural candidate. As pointed out in 
|4|], at the typical muon energy involved in un- 
derground analyses (E^ ~ 1 TeV), if very large 
energy transfer are considered, the cross section 
of such a process is non negligible with respect 
to the e + e~ pair production. The key point is 
that at high v—Ej/E^, the most important Feyn- 
mann diagram is rather insensitive to the ratio 
(me/m^) 2 . An analytic expression for the muon 
pair production cross section is given in j|,||. In 
order to test the hypothesis, such cross section 
has been included in the muon transport code 
PROPMUra]. Assuming a muon flux with energy 
spectrum E~ 3 ' 7 with a minimum muon energy 
£™ m = 1.2 TeV at the surface and considering 
the actual mountain profile, we generated a sam- 
ple of 10 muons corresponding to 3666 h of live 
time. About ~ 3.0 • 10 6 muons survived to the 
MACRO level, 5360 of which derived from muon 
pair production processes. The average distance 
of these muon pairs is (128 ± 1) cm and their aver- 
age energies are (657±14) GeV and (145±3) GcV 
for, the main muon and the secondary muon sam- 
ples respectively. We propagated the muons sur- 
viving to the MACRO level through GEANT sim- 
ulation and we applied to these events the same 
cuts specified in Section 2. Finally, the number 
of events was normalized to the live time of real 
data. 

In Table u\ we report the number of weighted 
muon pairs in the first bins of the experimen- 
tal and simulated decoherence distributions ex- 
pressed as percentage with respect to the popula- 
tion in the bin of the distribution maximum. The 
effect of standard cuts, of the C4 cut and of the 
subtraction of the muon pair production process 
are shown in order. In each case, we indicate in 
percentage the bin populations with respect to 
the peak of the distribution and the discrepancy 
with respect to the Monte Carlo predictions. 

In Fig. we compare the simulated decoher- 
ence curve with the experimental data corrected 
for the muon pair production effect. Despite the 
approximation introduced in our test, it seems 
that the proposed muon pair production process 
can account for the most of observed discrepancy 
in the low distance range. 



Table 1 

Number of weighted muon pairs in the first bins of the experimental and simulated decoherence dis- 
tributions. The discrepancy is the percentage difference between experimental and Monte Carlo values 
normalized to the distribution maximum (last column). 
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Figure 2. The low distance region of the experi- 
mental decoherence function before and after the 
subtraction of the secondary muon sample and 
comparison with the Monte Carlo simulation. 



3. Discussion and conclusion 



We measured the underground decoherence 
function using high energy muons (E^ > 
1.3 TeV) up to a maximum distance of about 
70 m. 

The capability to resolve very close muon 



tracks permitted the investigation of the be- 
haviour of the decoherence function for small sep- 
arations. Apart from the negligible contamina- 
tion of hadro-production by muons, we found 
strong hints of a contribution coming from muon 
pair production by muons. The inclusion of this 
process in the simulation reproduces, both in a 
qualitative and quantitative way, the experimen- 
tal data. 
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